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^ : ABSTRACT 
(N 

Physical properties of the quantum gravitational vacuum state are explored by solving a lattice 
^ ■ version of the Wheeler-DeWitt equation. The constraint of diffeomorphism invariance is strong 

, enough to uniquely determine the structure of the vacuum wave functional in the limit of infinitely 

fine triangulations of the three-sphere. In the large fluctuation regime the nature of the wave 
function solution is such that a physically acceptable ground state emerges, with a finite non- 
perturbative correlation length naturally cutting off any infrared divergences. The location of 
the critical point in Newton's constant Gc, separating the weak from the strong coupling phase, 
is obtained, and it is inferred from the structure of the wave functional that fluctuations in the 
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curvatures become unbounded at this point. Investigations of the vacuum wave functional further 
indicate that for weak enough couphng, G < Gc, a pathological ground state with no continuum 
limit appears, where configurations with small curvature have vanishingly small probability. One is 
then lead to the conclusion that the weak coupling, perturbative ground state of quantum gravity 
is non-perturbatively unstable, and that gravitational screening cannot be physically realized in 
the lattice theory. The results we find are in general agreement with the Euclidean lattice gravity 
results, and lend further support to the claim that the Lorentzian and Euclidean lattice formulations 
for gravity describe the same underlying non-perturbative physics. 



1 Introduction 



We have argued in previous work that the correct identification of the true ground state for quantum 
gravitation necessarily requires the introduction of a consistent nonperturbative cutoff, followed by 
the construction of the continuum limit in accordance with the methods of the renormalization 
group. To this day the only known way to introduce such a non-perturbative cutoff reliably in 
quantum field theory is via the lattice formulation. A wealth of results have been obtained over 
the years using the Euclidean lattice formulation, which allows the identification of the physical 
ground state and the accurate calculations of gravitational scaling dimensions, relevant for the scale 
dependence of Newton's constant in the universal scaling limit. 

In this work we will focus instead on the Hamiltonian approach to gravity, which assumes from 
the beginning a metric with Lorentzian signature. In order to obtain useful insights regarding 
the non-perturbative ground state, a Hamiltonian lattice formulation was written down based on 
the Wheeler-DeWitt equation, where the quantum gravity Hamiltonian is expressed in the metric- 
space representation. In [l] [2| a general discrete Wheeler-DeWitt equation was given for pure 
gravity, based on the simplicial lattice formulation of gravity developed by Regge and Wheeler. 
Here, we extend the work initiated in [U [2] to the physical case of 3 + 1 dimensions, and show 
how nonperturbative vacuum solutions to the lattice Wheeler-DeWitt equations can be obtained 
for arbitrary values of Newton's constant G. The procedure we follow is similar to what was 
done earlier in 2 + 1 dimensions. We solve the lattice equations exactly for several finite regular 
triangulations of the three-sphere, and then extend the result to an arbitrarily large number of 
tetrahedra. For large enough volumes the exact lattice wave functional is expected to depend 
on geometric quantities only, such as the total volumes and the total integrated curvature. The 
regularity condition on the solutions of the wave equation at small volumes is then shown to play 
an essential role in constraining the form of the vacuum wave functional. A key ingredient in 
the derivation of the results is the local diffeomorphism invariance of the Regge- Wheeler lattice 
formulation. 

From the structure of the resulting wave function a number of suggestive physical results can 
be obtained, the first one of which is that the non-perturbative correlation length is found to be 
finite for sufficiently large G. At the critical point G = Gc, which we determine exactly from the 
structure of the wave function, fiuctuations in the curvature become unbounded, thus signaling a 
divergence in the non-perturbative gravitational correlation length. Such a result can be viewed 



3 



as consistent with the existence of a non-trivial ultraviolet fixed point (or a phase transition in 
statistical field theory language) in G. The behavior of the theory in the vicinity of such a fixed 
point is then expected to determine, through standard renormalization group arguments, the scale 
dependence of the gravitational coupling in the vicinity of the ultraviolet fixed point. 

A short outline of the paper is as follows. In Sec. 2, as a general background to the rest of 
the paper, we briefly describe the formalism of canonical gravity. The continuum Wheeler-DeWitt 
equation and its invariance properties are introduced as well at this stage. Then we briefly outline 
the properties of the lattice Wheeler-DeWitt equation derived in our previous work, and later Sec. 3 
makes more explicit various quantities appearing in it. This section also emphasizes the important 
role of continuous lattice diffeomorphism invariance in the Regge framework, as it applies to the 
present case of 3 + 1-dimensional gravity. Sec. 4 focuses on the scaling properties of the lattice 
equations and various sensible choices for the lattice coupling constants, with the aim of giving 
eventually a more transparent form to the wave function results. Sec. 5 gives an outline of the 
general method of solution for the lattice equations, which are later discussed in some detail for 
a number of regular triangulations of the three-sphere. Then a general form of the wave function 
is given that covers all the previous discrete cases, and subsequently allows a study of the infinite 
volume limit. Sec. 6 discusses the issue of how to define an average volume and thus an average 
lattice spacing, an essential ingredient in the interpretation of the results given later. Sec. 7 
discusses modifications of the wave function solution obtained when the explicit curvature term in 
the Wheeler-DeWitt equation is added. Later a partial differential equation for the wave function 
is derived in the curvature and volume variables. General properties of the solution to this equation 
are discussed in Sec. 8. Sec. 9 contains a brief summary of the results obtained so far. 



2 Continuum and Discrete Wheeler-DeWitt Equation 



Since this paper involves the canonical quantization of gravity we begin here with a very brief 
summary of the classical canonical formalism [3j as derived by Arnowitt, Deser and Misner [^. 
While many of the results presented in this section are rather well known, it will be useful, in view 
of later applications, to recall the main results and formulas and provide suitable references for 
expressions used later in the paper. Here x* {i = 1,2, 3) will be coordinates on a three-dimensional 
manifold, and indices will be raised and lowered with gij{x) = 1,2,3), the three-metric on the 
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given spacelike hypersurface. As usual g^^ denotes the inverse of the matrix gij. Our conventions 
are such that the space-time metric has signature — h +- 1-, that ^i? is non- negative in a space-time 
containing normal matter, and that is positive in a 3-space of positive curvature. 

A transition from the classical to the quantum description of gravity is obtained by promoting 
the metric Qij, the conjugate momenta tt^^ , the Hamiltonian density H and the momentum density 
Hi to quantum operators, with gij and vr*-' satisfying canonical commutation relations. In particular, 
the classical constraints now select a physical vacuum state l^*), such that in the source- free case 

H\^) =0 H,\^) = , (1) 

and in the presence of sources more generally 

r|^') = fi\^) = , (2) 

where T and Tj now include matter contributions that should be added to H and Hi. As in ordinary 
nonrelativistic quantum mechanics, one can choose different representations for the canonically 
conjugate operators gij and vr*-'. In the functional metric representation one sets 

5,,(x) ^ 5y(x) 7r^J(x) ^ -ih-WnG-—^ . (3) 

In this picture quantum states become wave functionals of the three- metric gij(x), 

1^) ^ ^b.j(x)] . (4) 

The two quantum-constraint equations in Eq. ^ then become the Wheeler-DeWitt equation [5*, ^ 

and the momentum constraint listed below. In Eq. ([5]) Gij^ki is the inverse of the DeWitt superme- 
tric, given by 

Gij,kl = \9'^^'^ {gik9jl+ guQjk - QijOkl) ■ (6) 
The three-dimensional version of the DeWitt supermetric itself, G^^'^^x) is given by 

In the metric representation the diffeomorphism (or momentum) constraint reads 

|2i5.jVfe^ + Hf^ ^bii(x)] =0, (8) 
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where H'^ and Hf are possible matter contributions. In the fohowing, we shall set both of these to 
zero as we will focus here almost exclusively on the pure gravitational case. Then the last contraint 
represents the necessary and sufficient condition that the wave functional %p[g\ be an invariant under 
coordinate transformations [7J. 

We should also mention here that a number of basic issues need to be addressed before one can 
gain a full and consistent understanding of the dynamical content of the theory (see, for example, 
[9l dni O [131 El as a small set of representative references). These include possible problems of 
operator ordering, and the specification of a suitable Hilbert space, which entails at some point a 
choice for the norm of wave functionals, for example in the Schrodinger form 

|2 



ll^r = j dii[g] ngij] > (9) 

where djji[g\ is the appropriate (DeWitt) functional measure over the three-metric gij. In this work 
we will attempt to address those issues as they appear within the relevant calculations. 

The starting point for the following discussion will be the Wheeler-DeWitt equation for pure 
gravity in the absence of matter, Eq. ([5]), 

^•2 



(^>^OfG,,,.i.) ^^^J^^^_^ - - 2A)} *b.,(x)l = 



(10) 



and the diffeomorphism constraint of Eq. 



|2i<7^,(x)Vfc(x)^-^| ^'bii(x)] = 0. (11) 

Both of these equations express a constraint on the state |^') at every x, each of the form H{x) \ ^) = 
and Hi (x)|^') = 0. It is then natural to view Eq. (jlOp as made up of three terms, the first one 
identified as a kinetic term for the metric degrees of freedom, the second one involving — and 
thus seen as a potential energy contribution (of either sign, due to the nature of the 3-curvature 
■^i?), and finally the cosmological constant term proportional to + A ^ acting as a mass- like term. 
The kinetic term can be regarded as containing a Laplace-Beltrami-type operator acting on the 
6-dimensional Riemannian manifold of positive definite metrics gij, with Gij^ki acting as its con- 
travariant metric [7J . It was further shown in the quoted reference that the manifold in question has 
hyperbolic signature — + + with pure dilations of gij corresponding to timelike displacements 
within this manifold of metrics. 

Next we turn to the lattice theory. Here we will follow the procedure outlined in [Ij and 
discretize the continuum Wheeler-DeWitt equation directly, a procedure that makes sense in the 
lattice formulation, as these equations are still formulated in terms of geometric objects, for which 
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the Regge theory is very well suited. On a simplicial lattice [Ml [T71 [181 13 ISQl [ZD [12] (see for 
example [23], and references therein, for a more complete discussion of the Regge- Wheeler lattice 
formulation for gravity) one knows that deformations of the squared edge lengths are linearly related 
to deformations of the induced metric. In a given simplex a, take coordinates based at a vertex 0, 
with axes along the edges from 0. The other vertices are each at unit coordinate distance from 
(see Figure 1 as an example of this labeling for a tetrahedron). In terms of these coordinates, the 
metric within a simplex is given by 

9iji<T) = i (4 + /2._;2) . (12) 

Note that in the following discussion only edges and volumes along the spatial direction are in- 
volved. Then by a straightforward exercise of varying the squared volume of a given simplex a in 
d dimensions to quadratic order in the metric (in the continuum), or in the squared edge lengths 
belonging to that simplex (on the lattice), one is led to the key identification |241 [25] 

where the quantity G^^{P) is seen to be local, given that the sum over a only extends over those 
simplices which contain either the i or the j edge. In the formulation of [Ij it will be adequate to 
limit the sum in Eq. (1130 to a single tetrahedron, and define the quantity G^^ for that tetrahedron. 
Then, in schematic terms, the lattice Wheeler-DeWitt equation for pure gravity will have the form 

|-(167rG)2G.,(Z2)^|^ - [{I') - 2A] I ^l'] = , (14) 

with Gij{P) the inverse of the matrix G'''^{J?) given above. The range of the summation over i and 
j and the appropriate expression for the scalar curvature, in this equation, are discussed below and 
made explicit in Eq. ([T5|) . 

It is clear that Eqs. ([5]) or ([T^ express a constraint equation at each "point" in space. Indeed, 
the first term in Eq. ()14p contains derivatives with respect to edges i and j connected by a matrix 
element Gij which is nonzero only if i and j are close to each other, essentially nearest neighbor. One 
would therefore expect that the first term could be represented by just a sum of edge contributions, 
all from within one {d — l)-simplex a (a tetrahedron in three dimensions). The second term 
containing ^i?(/^) in Eq. (I14p is also local in the edge lengths: it only involves a handful of edge 
lengths, which enter into the definition of areas, volumes and angles around the point x. The latter 
is therefore described, through the deficit angle 5h, by a sum over contributions over all {d — 3)- 
dimensional hinges (edges in 3-1-1 dimensions) h attached to the simplex a. This then leads in three 
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dimensions to a more explicit form of Eq. ()14p 

< -(167rG)2 J2 G^ji<T)-gf^ - 2n^h Yl ^h^n + 2\vA -^[f] = . (15) 

i,jC(J * hCcT J 

Here 5h is the deficit angle at the hinge (edge) h, Ih the corresponding edge length, and = \/ g{(y) 
the volume of the simplex (tetrahedron in three spatial dimensions) labeled by a. Gij (a) is obtained 
either from Eq. (jl3p or from the lattice transcription of Eq. ([6]) 

Gij^kiicr) = \g~^''^{(T) [gik{(y)9ji{(y) + 9ii{(T)gjk{<y) - gij{<y)gki{<y)] , (16) 

with the induced metric gijic) within a simplex a given in Eq. p2p . The combinatorial factor rif^h 
ensures the correct normalization for the curvature term, since the latter has to give the lattice 
version of j ^ = 2 5hlh (in three spatial dimensions) when summed over all simplices a. The 
inverse of n^jh counts, therefore, the number of times the same hinge appears in various neighboring 
simplices and consequently depends on the specific choice of underlying lattice structure. The lattice 
Wheeler-DeWitt equation given in Eq. (jlSp was the main result of a previous paper pT] and was 
studied extensively in 2 + 1 dimensions in previous work |2]. 



3 Explicit Setup for the Lattice Wheeler-DeWitt Equation 

In the following we will now focus on a lattice made up of a large number of tetrahedra, with 
squared edge lengths considered as the fundamental degrees of freedom. For ease of notation, we 
define /qi = cl^ ^12 = b, Iq2 = c, /qs = d, /fa = e, /^s = /• the tetrahedron labeled as in Figure 1, 
we have 

gu = a, 522 = c, = d, (17) 

111 

512 = -(a + c - 6) , 513 = -{a + d - e), 523 = -^{c + d - f) , (18) 

and its volume V is given by 

V"^ = r^lafi-a- f + b + c + d + e) + bd{-b-d + a + c + e + f) + 

+ ce(— c — e + a + b + d+f) — abc — ade — bef — cdf] . (19) 
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The matrix G^^ is then given by 



24V 



( -2/ e + /-6 6 + /-e d + /-c c + /-d p \ 

e + f — b —2e b + e — f d + e — a q a + e — d 

b + f — e b -\- e — f —2b r b-\-c — a a + b — c 

d + f — c d + e — a r —2d c + d — f a + d — e 

c + f — d q b + c — a c + d — f —2c a + c — b 

\ p a + e — d a + b — c a + d — e a + c — b —2a / 



(20) 



where the three quantities p, q and r are defined as 
p = -2a -2f + b + c + d + e, q = -2c - 2e + a + b + d + f, r = -2b-2d + a + c + e + f. (21) 
To obtain Gij one can then either invert the above expression, or evaluate 

Gij^ki = i^—pigikgji + gugjk - gijgki), (22) 
g 

and later replace derivatives with respect to the metric by derivatives with respect to the squared 
edge lengths, as in ~ Wa ^ WB ^ We finds [Ij that the matrix representing the 

coefficients of the derivatives with respect to the squared edge lengths is the same as the inverse of 
C"^ , a result that provides a nontrivial confirmation of the correctness of the Lund-Regge result of 
Eq. Then in 3 + 1 dimensions, the discrete Wheeler-DeWitt equation is 

|-(16^G)2g,,^^ - 2n^hY.^6h + 2XV^^s] = 0, (23) 

where the sum is over hinges (edges) h in the tetrahedron. Note the mild nonlocality of the equation 
in that the curvature term, through the deficit angles, involves edge lengths from neighboring 
tetrahedra. In the continuum, the derivatives also give some mild nonlocality. Figure 2 gives a 
pictorial representation of lattices that can be used for numerical studies of quantum gravity in 
3+1 dimensions. 

In the following we will be concerned at some point with various discrete, but generally regular, 
triangulations of the three-sphere [2b\ [27] . These were already studied in some detail within the 
framework of the Regge theory in [20], where in particular the 5-cell 04, the 16-cell /?4 and the 
600-cell regular polytopes (as well as a few others) were considered in some detail. For a very early 
application of these regular triangulations to general relativity see [58]. 

We shall not dwell here on a well-known key aspect of the Regge- Wheeler theory, which is the 
presence of a continuous, local lattice diffeomorphism invariance, whose main aspects in regards to 
its relevance for the 3 + 1 formulation of gravity were already addressed in some detail in various 
works, both in the framework of the lattice weak field expansion |16[ |T], and beyond it |22[ 129] . 
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Figure 1: A tetrahedron with labels. 




Figure 2: A small section of a suitable spatial lattice for quantum gravity in 3 + 1 dimensions. 
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Here we will limit ourselves to some brief remarks on how this local invariance manifests itself in the 
3 + 1 formulation, and, in particular, in the case of the discrete triangulations of the sphere studied 
later on in this paper. In general, lattice diffeomorphisms in the Regge- Wheeler theory correspond 
to local deformations of the edge lengths about a vertex, which leave the local geometry physically 
unchanged, the latter being described by the values of local lattice operators corresponding to local 
volumes and curvatures [22t I29j . The case of flat space (curvature locally equal to zero) or 
near-flat space (curvature locally small) is obviously the simplest to analyze |29j: by moving the 
location of the vertices around on a smooth manifold one can find different assignments of edge 
lengths representing locally the same flat, or near- flat, geometry. Then it is easy to show that 
one obtains a d- A'o-parameter family of local transformations for the edge lengths, as expected for 
lattice diffeomorphisms. For the present case, the relevant lattice diffeomorphisms are the ones that 
apply to the three-dimensional, spatial theory. The reader is referred to [30] and, more recently, 
m for their explicit form within the framework of the lattice weak field expansion. 

With these observations in mind, we can now turn to a discussion of the solution method for the 
lattice Wheeler-DeWitt equation in 3 + 1 dimensions. One item that needs to be brought up at this 
point is the proper normalization of various terms (kinetic, cosmological and curvature) appearing 
in the lattice equation of Eqs. (jl5p and (j23p . For the lattice gravity action in d dimensions one has 
generally the following correspondence 

I'd'^x^ ^ J^K, (24) 

where Va is the volume of a simplex; in three dimensions it is simply the volume of a tetrahedron. 
The curvature term involves deficit angles in the discrete case, 

i / d^'x^R ^ Yl ^hSh , (25) 

h 

where 6fi is the deficit angle at the hinge h, and the associated "volume of a hinge" [15] . In four 
dimensions the latter is the area of a triangle (usually denoted by V^), whereas in three dimensions 
it is simply given by the length Ih of the edge labeled by h. In this work we will focus almost 
exclusively on the case of 3 + 1 dimensions; consequently the relevant formulas will be Eqs. (j24p 
and (I25p for dimension d = 3. 

The continuum Wheeler-DeWitt equation is local, as can be seen from Eq. (jlOp . One can 
integrate the Wheeler-DeWitt operator over all space and obtain 

|- {IGirGf j d^x A{g) + 2A j d-^x^- j d^x^R^ ^ = , (26) 
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with the super-Laplacian on metrics defined as 

A(<7) ^ G',,,h(x) ^ , , . . (27) 

In the discrete case one has one local Wheeler-DeWitt equation for each tetrahedron [see Eqs. ([T 
and (1151)]. which therefore takes the form 



[ hCa J 



(28) 



where A(/^) is the lattice version of the super-Laplacian, and we have set for convenience k = 2n^ h. 
As we shall see below, for a regular lattice of fixed coordination number, k is a constant and does 
not depend on the location on the lattice. In the above expression A(/^) is a discretized form of 
the covariant super-Laplacian, acting locally on the space of s = variables, 



with the matrix G^^ given explicitly in Eq. (j20p . Note that the curvature term involves six deficit 
angles 6h, associated with the six edges of a tetrahedron. Now, Eq. ([23]) applies to a single given 
tetrahedron (labeled here by cr), with one equation to be satisfied at each tetrahedron on the lattice. 
But one can also construct the total Hamiltonian by simply summing over all tetrahedra, which 
leads to 

I cr cr (J hCa ) 

Summing over all tetrahedra (cr) is different from summing over all hinges (h), and the above 
equation is equivalent to 



|-(167rG)2 ^A(;2) + 2A^y, - KqJ2lhSh \ ^ = 

I cr ah) 



(31) 



where q here is the lattice coordination number. The latter is determined by how the lattice is put 
together (which vertices are neighbors to each other, or, equivalently, by the so-called incidence 
matrix). Here, q is the number of neighboring simplexes that share a given hinge (edge). For a 
flat triangular lattice in 2d g = 6, whereas for the regular triangulations of we will be consid- 
ering below one has g = 3, 4, 5; for more general, irregular triangulations q might change locally 
throughout the lattice. For proper normalization in Eq. (I30p one requires the three-dimensional 
version of Eqs. ()24p and (I25p . which flxes the overall normalization of the curvature term 

K = 2n^h = -, (32) 
Q 

thus determining the relative weight of the local volume and curvature terms. 
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4 Choice of Coupling Constants 



We will find it convenient, in analogy to what is commonly done in the Euclidean lattice theory 
of gravity, to factor out an overall irrelevant length scale from the problem, and set the (unsealed) 
cosmological constant equal to one [20] . Indeed, recall that the Euclidean path integral statistical 
weight always contains a factor P{V) oc exp(— AqV), where V = J is the total volume on the 
lattice, and Aq is the unsealed cosmological constant. A simple global rescaling of the metric (or 
edge lengths) then allows one to entirely reabsorb this Aq into the local volume term. The choice 
Ao = 1 then trivially fixes this overall scale once and for all. Since Aq = 2A/167rG', one then has 
A = SttG in this system of units. In the following we will also find it convenient to introduce a 
scaled coupling A defined as 

A ^ (-^Y . (33) 

Then for Aq = 1 (in units of the UV cutoff or, equivalently, in units of the fundamental lattice 
spacing) one has A = l/647rG. 

Two further notational simplifications will be useful in the following. The first one is introduced 
in order to avoid lots of factors of IGvr in many of the formulas. So from now on we shall write G 
as a shorthand for IGvr G, 

I&ttG^G . (34) 

In this new notation one has X = G/2 and A = 1/4G. The above notational choices then lead to a 
more streamlined representation of the Wheeler-DeWitt equation, namely 

|-A + 1^5- ^^^i?} ^ = 0. (35) 

Note that we have arranged things so that now the kinetic term (the term involving the Laplacian) 
has a unit coefficient. Then in the extreme strong coupling limit (G — t- oo) the kinetic term is the 
dominant one, followed by the volume (cosmological constant) term (using the facts about A given 
above) and, finally, by the curvature term. Consequently, at least in a first approximation, the 
curvature R term can be neglected compared to the other two terms, in this limit. 

A second notational choice will later be dictated by the structure of the wave function solutions, 
which often involve numerous factors of VG. It will therefore be useful to define a new coupling g 
as 

g = VG, (36) 
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so that A = 4/(7^ (the latter g should not be confused with the square root of the determinant of 
the metric). 

5 Outline of the General Method of Solution 

The previous discussion shows that in the strong coupling limit (large G) one can, at least in 
a first approximation, neglect the curvature term, which will then be included at a later stage. 
This simplifies the problem considerably, as it is the curvature term that introduces complicated 
interactions between neighboring simplices. 

Here the general procedure for finding a solution will be rather similar to what was done in 2 + 1 
dimensions, as the formal issues in obtaining a solution are not dramatically different. First an 
exact solution is found for equilateral edge lengths s. Later this solution is extended to determine 
whether it is consistent to higher order in the weak field expansion, where one writes for the squared 
edge lengths the expansion 

= s{l + ehij) , (37) 

with e a small expansion parameter. The resulting solution for the wave function can then be 
obtained as a suitable power series in the h variables, combined with the standard Frobenius 
method, appropriate for the study of quantum mechanical wave equations for suitably well-behaved 
potentials. In this method one first determines the correct asymptotic behavior of the solution for 
small and large arguments, and later constructs a full solution by writing the remainder as a power 
series or polynomial in the relevant variable. While this last method is rather time consuming, we 
have found nevertheless that in some cases (such as the single triangle in 2 + 1 dimensions and the 
single tetrahedron in 3 + 1 dimensions, described in [1] and also below), one is lucky enough to find 
immediately an exact solution, without having to rely in any way on the weak field expansion. 

More importantly, in [2\ it was found that already in 2 + 1 dimensions this rather laborious 
weak field expansion of the solution is not really necessary, for the following reason. Diffeomorphism 
invariance (on the lattice and in the continuum) of the theory severely restricts the form of the 
Wheeler-DeWitt wave function to a function of invariants only, such as total three-volumes and 
curvatures, or powers thereof. In other words, the wave function is found to be a function of 
invariants such as / d'^x^/g or J d'^x^ etc. (these will be listed in more detail below for the 
specific case of 3 + 1 dimensions, where one has d = 3 in the above expressions). 
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For concreteness and computational expedience, in the following we will look at a variety of 
three-dimensional simplicial lattices, including regular triangulations of the three-sphere con- 
structed as convex 4-polytopes, the latter describing closed and connected figures composed of lower 
dimensional simplices. Here these will include the 5-cell 4-simplex or hypertetrahedron (Schlafli 
symbol {3,3,3}) with 5 vertices, 10 edges and 5 tetrahedra; the 16-cell hyperoctahedron (Schlafli 
symbol {3,3,4}) with 8 vertices, 24 edges and 16 tetrahedra; and the 600-cell hypericosahedron 
(Schlafli symbol {3, 3, 5}) with 120 vertices, 720 edges and 600 tetrahedra l26l[27]. Note that the Eu- 
ler characteristic for all 4-polytopes that are topological 3-spheres is zero, x = A'^q— -^i+-^2— -^3 = 0, 
where A'^^ is the number of simplices of dimension d. We also note here that there are no known reg- 
ular equilateral triangulations of the flat 3-torus in three dimensions, although very useful slightly 
irregular (but periodic) triangulations are easily constructed by subdividing cubes on a square 
lattice into tetrahedra [30J. 

In the following we will also recognize that there are natural sets of variables for displaying the 
results. One of them is the scaled total volume x, defined as 

_ 4V2X ^ 4V2X 

a 

Later on we will be interested in making contact with continuum manifolds, by taking the infinite 
volume (or thermodynamic) limit, defined in the usual way as 

A^o- — ^ oo , 
Vtot oo , 

— ^ — >• const. , (39) 

with N(j = N3 here the total number of tetrahedra. It should be clear that this last ratio can be 
used to define a fundamental lattice spacing ao, for example via Vtot/^a = Va = ciq/gV^- 

The full solution of the quantum mechanical problem will, in general, require that the wave 
functions be properly normalized, as in Eq. ([9]). This will introduce at some stage wave function 
normalization factors M, which will later be fixed by the standard rules of quantum mechanics. 
If the wave function were to depend on the total volume Vtot only (which is the case in 2 + 1 
dimensions, but not in 3 + 1), then the relevant requirement would simply be 

/poo 
df^[9]-\mj]\^= dVtofVt'^t-\'^iytot)\^= 1, (40) 

where d^[g] is the appropriate functional measure over the three- metric gij, and m a positive real 
number representing the correct entropy weighting. But, not unexpectedly, in 3 + 1 dimensions 
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the total curvature also plays a role, so the above can only be regarded as roughly correct in the 
strong coupling limit (large G), where the curvature contribution to the Wheeler-DeWitt equation 
can safely be neglected. As in nonrelativistic quantum mechanics, the normalization condition in 
Eqs. Q and (j40p plays a crucial role in selecting out of the two solutions the one that is regular, 
and therefore satisfies the required wave function normalizability condition. 

To proceed further, it will be necessary to discuss each lattice separately in some detail. For 
each lattice geometry, we will break down the presentation into two separate discussions. The first 
part will deal with the case of no explicit curvature term in the Wheeler-DeWitt equation. Each 
regular triangulation of the three-sphere will be first analyzed separately, and subjected to the 
required regularity conditions. Here a solution is first obtained in the equilateral case, and later 
promoted on the basis of lattice diffeomorphism invariance to the case of arbitrary edge lengths, 
as was done in Later a single general solution will be written down, involving the parameter 
q, which covers all previous triangulation cases, and thereby allows a first study of the infinite 
volume limit. The second part deals with the extension of the previous solutions to the case when 
the curvature term in the Wheeler-DeWitt equation is included. This case is more challenging to 
treat analytically, and the only results we have obtained so far deal with the large volume limit, for 
which the solution is nevertheless expected to be exact (as was the case in 2 + 1 dimensions [2]). 

5.1 Nature of Solutions in 3+1 Dimensions 

In this work we will be concerned with the solution of the Wheeler-DeWitt equation for discrete 
triangulations of the three-sphere S"^. In general, for an arbitrary triangulation of a smooth closed 
manifold in three dimensions, one can write down the Euler equation 

A^o - iVi + iV2 - iVs = (41) 

and the Dehn-Sommerville relation 

N2 = 2N3. (42) 

The latter follows from the fact that each triangle is shared by two tetrahedra and each tetrahedron 
has four triangles, thus 2N2 = 4 A'^s. In addition, for the regular triangulations of the three-sphere 
we will be considering here, one has the additional identity 

A^i = - iV3 , (43) 

q 
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where q is the local coordination number, defined as the number of tetrahedra meeting at an edge. 
For the three regular triangulations of the three-sphere we will look at one has g = 3, 4, 5. The 
above relations then allow us to relate the number of sites (iVg) to the number of tetrahedra {N^)^ 

No = iVg (^^ - 1^ . (44) 

It will also turn out to be convenient to collect here a number of useful definitions, results and 
identities that apply to the regular triangulations of the three-sphere, valid strictly when all edge 
lengths take on the same identical value I = y/s. For the total volume 

Vtot= Y.^- ^ j (45) 

one has 

^3/2 



whereas the total curvature 



Vtot = iVs K = ^ iVs , (46) 



Rtot = 2Y,^hlh ^ I d^x^R (47) 



h 

is given by 

Rtot = ^ [2 7r-gcos-i(|)] iVs. (48) 
The latter relationship can be inverted to give the parameter g as a function of the curvature 

10 ^^^^\, (49) 



and its inverse as 



Rtot + ^ iV3 



-1 _ -1 , Rtot 



SO that this last quantity is just linear in Rtot- A very special value for q corresponds to the choice 
q = qo for which Rtot = 0. For this case one has 

qo = = 5.1043 . (51) 

cos^^(^) 

Then, summarizing all the previous discussions, the discretized Wheeler-DeWitt equation one wants 
to solve here is 

Gii(^)^|^ - ^ E ^h5h + 2\vA = , (52) 

i,jGa- * i hCcr J 

with parameter k given by 

K = ^ . (53) 

If the reader is not interested in the details of the solution for each individual lattice, then (s)he 
can skip the following sections and move on directly to Sec. (15. 6p . 
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5.2 1-Cell Complex (Single Tetrahedron) 

As a first case we consider here the quantum-mechanical problem of a single tetrahedron. One 
has A'o = 4, A'^i = 6, = 4, = 1 and q = 1 [note that these do not satisfy the Euler and 
Dehn-Sommerville relations; only the relation between A'^i , N^, and q, Eq. (|44p . is satisfied for a 
single tetrahedron]. The single tetrahedron problem is relevant for the strong coupling (large G) 
limit. In this limit one can neglect the curvature term, which couples different tetrahedra to each 
other, and one is left with the local degrees of freedom, involving a single tetrahedron. 

The Wheeler-DeWitt equation for a single tetrahedron with a constant curvature density term 
R reads 

I - (167rG)2 Gij-^^ + (2A - i?) y| vl/[s] = , (54) 

where now the squared edge lengths si . . . sq are all part of the same tetrahedron, and Gij is given 
by a rather complicated, but explicit, 6x6 matrix given earlier. 

As in the 2+1 case previously discussed in [2], here too it is found that, when acting on functions 
of the tetrahedron volume, the Laplacian term still returns some other function of the volume only, 
which makes it possible to readily obtain a full solution for the wave function. In terms of the 
volume of the tetrahedron one has the equivalent equation for ^[s] = ^(Va) (note that we have 
now replaced for notational convenience 167rG — )• G) 

7 S2 A 

i>" (T4) + - V' {V.) + (K) = , (55) 

with primes indicating derivatives with respect to 14-. From now on we will set the constant 
curvature density R=0. If one introduces the dimensionless (scaled volume) variable 

4\/2A 

X = —Q— ^tot , (56) 

where Vtot = is the volume of the tetrahedron, then the differential equation for a single tetra- 
hedron becomes simply 

i;" (x) + - iIj' (x) + 4^ (x) = 0. (57) 

X 

Solutions to Eqs. ([5^ or ([FT]) are Bessel functions Jm or Ym with m = 3, 



i^RiVtot) = const. J3 I ^^Vtot J /Vtlt , (58) 



or 



MVtot) = const. I ^^Vtot ) /V-it . (59) 
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Only Jrnix) is regular as x — >• 0, Jm{x) ~ V{m + 1) In terms of the variable x the regular 

solution is therefore 

and the only physically acceptable wave function is 



V' {Vtot) oc cc ^ ^ 3 , (60) 



T /4v12At/ 
-'3 —(7— l^tot 

f (a, 6, . . . /) = ^(Ftoi) = , (61) 

not 

with normalization constant 

^ " 1024 23/4 [^) ■ ^^^^ 
The latter is obtained from the wave function normalization requirement 

) 

dVtot I ^{Vtot) I' = 1 . (63) 



Note that the solution given in Eq. (j60p is exact, and a function of the volume of the tetrahedron 
only; its only dependence on the values of the edge lengths of the tetrahedron [or, equivalently, on 
the metric, see Eq. (jl2|) ] is through the total volume. 

One can compute the average volume of the single tetrahedron, and it is given by 



( Vtot > 



r dVtot ■ Vtot ■ I "^{Vtot) ? = ^^^^^""^ = 0.2643 . (64) 
Jo ' ^ 262144 V2A ^/\ ^ ' 



This last result allows us to define an average lattice spacing, by comparing it to the value for an 
equilateral tetrahedron for which Vtot = (l/6\/2) Oq- One obtains 

/ (-r X 1/3 

ao = 1.3089 ^-^J . (65) 

In terms of the parameter A defined in Eq. ([33]) one has \/\/G = \f2\. With the notation of 
Eq. p6|) one has as well G/\fX = yj2G = \f2g. Then for a single tetrahedron one has ( Vtot ) = 
(V^) = 0.37385. 

The single tetrahedron problem is clearly quite relevant for the limit of strong gravitational 
coupling, 1 /G — )• 0. In this limit lattice quantum gravity has a finite correlation length, comparable 
to one lattice spacing, 

^ ~ ao . (66) 

This last result is seen here simply as a reflection of the fact that for large G the edge lengths, and 
therefore the metric, fluctuate more or less independently in different spatial regions, due to the 
absence of the curvature term in the Wheeler-DeWitt equation. This is of course true also in the 
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Euclidean lattice theory, in the same limit [20] . It is the inclusion of the curvature term that later 
leads to a coupling between fluctuations in different spatial regions, an essential ingredient of the 
full theory. 

5.3 5-cell Complex (Configuration of 5 Tetrahedra) 

The first regular triangulation of S"^ we will consider is the 5-cell complex, sometimes referred to as 
the hypertetrahedron. Here one has Nq = 5, Ni = 10, = 10, N3 = 5 and q = 3, since there are 
three tetrahedra meeting on each edge. Then for the parameter k appearing in Eq. (j52p one has 



First we will consider the case of no curvature term in the lattice Wheeler-DeWitt equation of 
Eq. (j52p . The curvature term will be re- introduced at a later stage [see Sec. d?])], as its presence 
considerably complicates the solution of the lattice equations. 

Solving the lattice equations directly (by brute force, one might say) in terms of the edge length 
variables is a rather difficult task, since many edge lengths are involved, increasingly more so for 
finer triangulations. Nevertheless it can be done, to some extent, in 2 + 1 dimensions [2], and 
possibly even in 3 + 1 dimensions, analytically for some special cases, or numerically for more 
general cases. To obtain a full solution to the lattice equations we rely here instead on a simpler 
procedure, already employed successfully (and checked explicitly) in 2 + 1 dimensions. First, an 
exact wave function solution to the lattice Wheeler-DeWitt equations is obtained for the equilateral 
case, where all edges in the simplicial complex are assumed to have the same length. Then, in the 
next step, the diffeomorphism invariance of the simplicial lattice theory is used to promote the 
previously obtained expression for the wave function to its unique general coordinate invariant 
form, involving various geometric volume and curvature terms. It is a non-trivial consequence of 
the invariance properties of the theory that such an invariant expression can be obtained, without 
any further ambiguity. In a number of instances such a procedure can be checked explicitly and 
systematically within the framework of the weak field expansion, and used to show that the form of 
the relevant wave function solution is indeed, as expected, strongly constrained by diffeomorphism 
invariance [2j. 

In the case of the 5-cell complex, and for now without an explicit curvature term in the Wheeler- 
DeWitt equation, one obtains the following differential equation 



2 



(67) 




V'" (Vtot) + 



95 



V'' (Vtot) + 



32 A 



V' (Vtot) = 



(68) 



914, 



tot 



9G2 
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for a wave function that, for now, depends only on the total volume, ij: = ip {Vtot)- To obtain this 
result, it is assumed at first that the simplicial complex is built out of equilateral tetrahedra; in 
accordance with the previous discussion, this constraint will be removed below. In terms of the 
dimensionless variable x defined as 



4V2A , , 

x^^^^toi (69) 



one has the equivalent form for Eq. (|68 



^l,"{x) + ^ij'{x) + ^l,{x) = (70) 
9 X 



This last equation can then be solved immediately, and the solution is 



J43 ix) J43 { ^nf/" Vtot 

i'iVtot) oc « ^ 43 , (71) 

Vtot' 

up to an overall wave function normalization constant. As in the previously discussed tetrahedron 
case, and also as in 2 + 1 dimensions, one discards the Bessel function of the second kind (Y) 
solution, since it is singular at the origin. 

5.4 16-cell Complex (Configuration of 16 Tetrahedra) 

The next regular triangulation of S"^ we will consider is the 16-cell complex, sometimes referred to 
as the hyperoctahedron. One has in this case Nq = 8, Ni = 24, = 32, N3 = 16 and (7 = 4, since 
there are four tetrahedra meeting on each edge. For the parameter k in Eq. (|52p one has 

^ = 1- (72) 

In the case of the 16-cell complex (again for now without an explicit curvature term in the Wheeler- 
DeWitt equation) one obtains the following differential equation 

47 2 A 

^" iVtot) + TTf^ V'' (Vtot) + 7^ ^ (Vtot) = (73) 

^ Vtot 

for a wave function that depends only on the total volume, tp = tp iytot)- In terms of the variable 

X = ^ Vtot (74) 



one has an equivalent form for Eq. (|73 



47 

i^" {x) + —i^' {x) + ilj{x) = Q . (75) 
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The correct wave function solution is now 

J45 (x) J45 i^^Vtot 
^ iVtot) oc oc ' ^ ^ , (76) 

up to an overall wave function normalization constant. Again, we discarded the Bessel function of 
the second kind {Y) solution, since it is singular at the origin. 

5.5 600-cell Complex (Configuration of 600 Tetrahedra) 

The last, and densest, regular triangulation of S"^ we will consider here is the 600-cell complex, 
often called the hypericosahedron. For this lattice one has Nq = 120, A'^i = 720, = 1200, 
= 600 and q = 5, since there are now five tetrahedra meeting at each edge. For the parameter 
n in Eq. ()52p one has 

^ = 1 (77) 

For this 600-cell complex (again for now without an explicit curvature term in the Wheeler-DeWitt 
equation) one obtains the following differential equation 

^" (Vtot) + ^ V'' {Vtot) + (Vtot) = (78) 

for a wave function that depends only on the total volume, ip = tp (Vtot)- In terms of the variable 

4\/2A 



^ - 5^ Vtot (79) 
one has an equivalent form for Eq. ()78p 

672 

^"{x) + V'(^) + V'(^) = . (80) 

X 

Then the solution of the Wheeler DeWitt equation without a curvature term is 

JeTi (x) Jen (^-^^ Vtot 
i> {Vtot) oc « ^ ^ 671 , (81) 

Vtot ' 

again up to an overall wave function normalization constant. As in previous cases, we discard the 
Bessel function of the second kind (Y) solution, since it is singular at the origin. 

5.6 Summary and General Case for Zero Curvature 

In this section we summarize and extend the previous results for the wave functions, obtained so 
far for the three separate cases of the 5-cell, 16-cell and 600-cell triangulation of the three-sphere 
. The single tetrahedron case is somewhat special (it cannot contain a curvature term), and will 
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be left aside for the moment. Also, all the previous results so far apply to the case of no explicit 
curvature term in the Wheeler-DeWitt equation of Eq. (j52p : the inclusion of the curvature term 
will be discussed later. Consequently the following discussion still focuses on the strong coupling 
limit, G —7- oo. 

For the following discussion the relevant Wheeler-DeWitt equation is the one in Eq. (j52p . 



-G^ Gii(a)^^ - ^ E ^^^^ + 2Ay. i = , (82) 

iJCu * J hGer J 

which depends on the parameter 

K = - , (83) 

where q represents the number of tetrahedra meeting at an edge. The above equation is quite 

general and not approximate in any way. Nevertheless it depends on the local lattice coordination 

number q (how the edges are connected to each other, or, in other words, on the incidence matrix). 

Now, all previous differential equations for the wave function as a function of the total volume 

Vtot [Eqs. ([68]) . ([73]) and ([78]) ] can be summarized as a single equation 

, , (11 + 9g) TVs , , , 32 A , , 
2 q^ Vtot G^ 

Equivalently, in terms of the scaled volume variable defined as 

4\/2A 



Vtot , (85) 



qG 

one can summarize the results of Eqs. ([70]) . ([75]) and ([80]) through the single equation 



^P"ix) + ^^^^^^—^'ix) + ijix) = 0. (86) 
2 X 



It will be convenient here to define the (Bessel function) index n as 



n + 9q 

^ - A 2 


1 

" 2 ' 




(87) 


so that for the 5-cell, 16-cell and 600-cell one has 








95 , 

2n + l = y (g = 3. 


iVs = 


5), 




= y (« = 4, 


iV3 = 


16), 




= 672 {q = 5, 


iV3 = 


600) , 


(88) 


respectively, and in the general case 








2 






(89) 
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thus reproducing n = 43/9, 45/4 and 671/2, respectively, in the three cases. Then Eq. ([86]) is just 

2n + 1 



^" {x) + 



ip' (x) + ip (x) 







(90) 



Consequently the wave function solutions are 

ip oc 



oc 



J„ (x) -^ri [ Vtot 



4V2A T/ 



(91) 



up to an overall wave function normalization constant, thus summarizing all the results so far for the 
individual regular triangulations [Eqs. (fTTI) . (i76|) and (fSTI) ]. A more explicit, but less transparent, 
form for the wave function solution is 



tot 



l_JV3{ll+9<j) 



1 , iV3(ll+9q) 
o H 7—7. 



4g 



4V2A 
q G 



Vt 



tot 



(92) 



with N an overall wave function normalization constant. Its large volume behavior is completely 
determined by the asymptotic expansion of the Bessel J function. 



il^^x) 



Jn(x) 



~ X 
x — 00 



— sm(x+- ) + 

irx 



1 



(93) 



4 2 

It is also easy to see that the argument of the Bessel function solution J in Eqs. (j9ip and (j92p has 
the following expansion for large volumes 



4 V2A 



X 



Vtot + 



R. 



tot 1 



(94) 



qoG ' 36^/2^ G 

with ao (cq = 6-v/2 V^/A^s) representing here the average lattice spacing. Thus the second correction 
is of order iV/N^)'^/'^ Rtot- Note that nothing particularly interesting is happening in the structure 
of the wave function so far. Similarly, the index n of the Bessel function solution in Eqs. (19ip and 
(|92p has the following expansion for large volumes and small curvatures, 

(22 + 9(?o) 



(11 + 9 go) 1 

n = r5 JS^ 

4a^ 2 



Rtot + R' 



(95) 



96 vr Qo ao 

with oo again defined as above. Note here that the second correction is of order {N^/Vy^^ Rtot- It 
follows that the asymptotic behavior for the exponent of the fundamental wave function solutions 
for large volume and small curvature is given by 



± i 



4V2A 
1^ 



Vtot + 



36v/2 7r G 
' ll + 9go 
Wo 



Rtot + O {R'] 

22 + 9go 



N3 + 



96 vr go oo 



Rtot + O (i?2) 



InVt 



tot 



(96) 
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Let us make here some additional comments. One might wonder what concrete lattices corre- 
spond to values of n greater that 671/2, which is after all the highest value attained for a regular 
triangulation of the three-sphere, namely the 600-cell complex. For each of the three regular tri- 



of tetrahedra meeting at an edge (the local coordination number). In the three cases examined 
previously q was of course an integer between three and five; in two dimensions it is possible to 
have one more integer value of q corresponding to the regularly triangulated torus, but this is not 
possible here. In any case, one always has for a given triangulation of the three-sphere the Euler 
relation A'^o — Ni + N2 — = 0. The interpretation of other, even noninteger, values of q is then 
clear. Additional triangulations of the three-sphere can be constructed by considering irregular 
triangulations, where the parameter q is now seen as an average coordination number. Of course 
the simplest example is what could be described as a semiregular lattice, with Na edges with coor- 
dination number qa and edges with coordination number q^, such that Na + Nh = Ni. Various 
irregular and random lattices were considered in detail some time ago in [19], and we refer the 
reader to this work for a clear exposition of the properties of these kind of lattices. In the following 
we will assume that such constructions are generally possible, so that even non-integer values of q 
are meaningful and are worth considering. 

6 Average Volume and Average Lattice Spacing 

At this stage it will be useful to examine the question of what values are allowed for the average 
volume. The latter will be needed later on to give meaning to the notion of an average lattice 
spacing. In general the average volume is defined as 



where dfj,[g] is the appropriate (DeWitt) functional measure over three-metrics gij. 

Now consider the wave function obtained given in Eq. (|9ip . with n defined in Eq. (j87|) . This 
wave function is relevant for the strong coupling limit, where the explicit curvature term in the 
Wheeler-DeWitt equation can be neglected. In this limit one can then compute the average total 
volume 



angulations of with A^'o sites one has for the number of edges A''i = q^Nq, for the number 
of triangles A'^2 = ^^^0 and for the number of tetrahedra = -^^Nq, where q is the number 




(97) 




(98) 
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One then obtains immediately for the average volume of a tetrahedron 

r (n) N3 V A 

If the whole lattice is just a single tetrahedron, then one has n = 3, and therefore 

, , 31185 vrG G , , 

from which one can define an average lattice spacing ao via (y^) = a|]/6 ^2. For large iVg one has 



^ U^^a, (101) 

But in general one cannot assume a trivial entropy factor from the functional measure, and one 
should evaluate instead 

with some power m = cq and cq a real positive constant. One then obtains for the average 
volume of a single tetrahedron 

= {Vtot) = Vco[n + qo{9-Coqo)] , (103) 

which is finite as A^3 — )• 00. Note that in order for the above expression to make sense one requires 
Co < (ll + 9go)/go — 2.185. If the exponent in the entropy factor is too large, the integrals diverge. 
One then finds that the corresponding lattice spacing is given by 

3G 



al = Vco[ll + (7o(9-co<Zo)] -y= ■ (104) 

4 vA 

The lesson learned from this exercise is that in gravity the lattice spacing ao (the fundamental 
length scale, or the ultraviolet cutoff if one wishes) is itself dynamical, and thus set by the bare 
values of G and A. In a system of units for which Ao = 1 one then has ao ~ g^^^- Either way, 
the choice for ao has no immediate direct physical meaning, and has to be viewed instead in the 
context of a subsequent consistent renormalization procedure. In the following it will be safe to 
assume, based on the results of Eqs. ()65p and (|104p that 

«o = /'^, (105) 

in units of the UV cutoff, where / is a numerical constant of order one (for concreteness, in the 
single tetrahedron case one has / 1.3089). 
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7 Large Volume Solution for Nonzero Curvature 

The next task in line is to determine the form of the wave function when the curvature term in the 
Wheeler-DeWitt equation of Eq. (j52p is not zero. In particular we will be interested in the changes 
to the wave function given in in Eqs. (j9ip and (j92p . with argument x in Eq. (|94p and parameter n 
in Eq. (f95|) . We define here the total integrated curvature Rtot as in Eq. ()T7|) . which is of course 
different from the local curvature appearing in the lattice Wheeler DeWitt equation of Eq. (j52p . 



R^ = Y,^hlh- (106) 

hCa- 

In order to establish the structure of the solutions for large volumes Vtot we will assume, based in 
part on the results of the previous sections, and on the analogous calculation in 2 + 1 dimensions 
[2], that the fundamental wave function solutions for large volumes have the form 



exp|±i^a j j d^x^R + j J d^Xy/^R^ + 5 j d^x^ Rf.^R'"' + 



■ (107) 



Note here that the structure of the above expression, and the nature of the terms that enter into it, 
are basically dictated by the requirement of diffeomorphism invariance as it applies to the argument 
of the wave functional. Apart from the cosmological term, allowed terms are all the ones that can 
be constructed from the Riemann tensor and its covariant derivatives, for a a fixed topology of 3- 
space. Clearly, at large distances (infrared limit) the most important terms will be the Einstein and 
cosmological terms, with coefficients /3 and a, respectively. In three dimensions the Riemann and 
Ricci tensor have the same number of algebraically independent components (6), and are related 
to each other by 

R'\. = ^'"^^x^P {R". - h^'.) ■ (108) 

The Weyl tensor vanishes identically, and one has 

Rf^uXaR^"^'' - ^R^.uR^" -3R^ = C^.XaC^"'^'' = . (109) 

As a consequence, there is in fact only one local curvature squared term one can write down in 
three spatial dimensions. Nevertheless, higher derivative terms will only become relevant at very 
short distances, comparable or smaller than the Planck length ^fG] in the scaling limit it is expectd 
that these can be safely neglected. 

When expressed in lattice language, the above form translates to an ansatz of the form 

exp{±i(coytot + ciRZ)} , (110) 
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with m assumed to be an integer. In addition, from the studies of lattice gravity 2 + 1 dimensions 
one expects a InVtot term as well in the argument of the exponential [2\. This suggests a slightly 
more general ansatz , 

exp {± i {coVtot + ciRZ) + C2lnVtot + C3 In Rtot} ■ (HI) 



The next step is to insert the above expression into the lattice Wheeler-DeWitt equation Eq. (j52p 
and determine the values of the five constants cq . . . C3, m. This can be done consistently just to 
leading order in the weak field expansion of Eq. (j37p . which is entirely adequate here, as it will 
provide enough information to uniquely determine the coefficients. Here we will just give the result 
of this exercise. For the 5-cell complex (q = 3) one obtains 

whereas for 16-cell complex (g = 4)one finds 

~ exp <^ ±d Vtot - Rtot j - -J InVtot \ , (113) 

and finally for 600-cell complex (q = 5) 

~ exp |± i 1^ ^^^^ Vtot - Rtot ^ - 336 In Vtot | • (114) 

These expressions allow us again to identify the answer for general q as 

/ / 4^/2A 3^/2 ^ \ (ll + 9g)jV3 , \ 

ip exp <±i [ — Vtot Rtot niVtot } ■ (115 

[ \ qGVX J J 

Note that in deriving the above results we considered the large volume limit V ^ 00, treating the 
number of tetrahedra as a fixed parameter. Then from the previous expression we can now read 
off the values for the various coefficients, namely 

4.V2X 

Co = 



ci 



C2 



qG 

' qQy/X 
(ll + 9(7)iV3 



Aq' 

C3 = (116) 



with the only possible value m = 1. 
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In order to make contact with the strong couphng result for the wave function derived in the 
previous sections [Eqs. ([92]) . (p5|) and ([96j) ]. one needs to again expand the above answer for 

smah curvatures. One obtains for the exponent of the wave function the following expression 



11 + 9(70 



iVs + 11^^ Rtot + O {R') ] In Vtot , (117) 
96 vr go ao J 



with oq again representing the average lattice spacing, Oq = 6\/^V/N3. This finally determines 
uniquely the coefficients a and f3 appearing in Eq. (|107p . 

4 V2X 



a 



* a 



^ = ^.^-A 1 . ,118) 

The most important result so far is the appearance of two contributions of opposite sign in /3, 
signaling the appearance of a critical value for G where f3 vanishes. 

This critical point is located at Ac = 108 \/27r/go ™ ^ system of units where A = G/2 

[see Sec. (g])], 

G. = ?15^4. (119) 

But since the average lattice spacing ao is itself a function of G and A [see Eqs. ([65|) . (jl04p and 
(jlOSp ] one obtains in the same system of units 

36 23/8 31/4 ^3/4 

Gc = ^3,, 3/4 - 28.512 , (120) 

using the value of / for the single tetrahedron, or equivalently Qc — 5.340, a rather large value. 
Nevertheless we should keep in mind that in this paper we are also using a system of units where 
we set IGvrG — )• G. So, in a conventional system of units, one has the more reasonable result 
Gc ~ 0.567 in units of the fundamental UV cutoff. Evidence for a phase transition in lattice 
gravity in 3 + 1 dimensions was also seen earlier from the application of the variational method, 
using Jastrow-Slater correlated product variational wavefunctions [1]. Q Note that the results of 
Eqs. pi7p and piSp imply a dependence of the fundamental wave function on the curvature, of 
the type 

iP{R) ~ e^^-f^*°*/^o , (121) 



^ One can compare the above value for Gc obtained in the Lorentzian 3 + 1 theory with the corresponding value 
in the Euclidean four-dimensional theory. There one finds Gc ~ 0.624 [32|. The two values are not expected to be 
the same in the two formulations, due to the different nature of the cutoffs. In particular, in the lattice Hamiltonian 
formulation the continuum limit has already been taken in the time direction. Nevertheless, it is encouraging that 
they are quite comparable in magnitude. 
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with Rq a characteristic scale for the total, integrated curvature. Thus Rq ~ l/{g — g^.) with Gc, 
and therefore Qc = ^/Gc, given in Eq. (|119p . Therefore at the critical point fluctuations in the 
curvature become unbounded, just as is the case for the fluctuations in a scalar field when the 
renormalized mass approaches zero. 

At this stage one can start to compare with the results obtained previously without the explicit 
curvature term in the Wheeler-DeWitt equation, Eqs. (j94p and (j95p . The main change is that here 
one would be led to identify 

^^\„J^^^_±.^]n,^, ,122) 



qoG \^36\/2 7r G qo GV2X^ 

so that the Bessel function argument x [see Eq. (I94p ] now contains a new contribution, of opposite 
sign, in the curvature term. Its origin can be traced back to the new curvature contribution ci 
in Eq. (jll6p . which in turn arises because of the explicit curvature term now present in the full 
Wheeler-DeWitt equation. On the other hand, as is already clear from the result for C2 in Eq. (I116p . 
the index n of the Bessel function solution in Eqs. (j9ip and ()92p is left unchanged, 

11 + 9^ 1 22 + 9go p ^n(T^2 \ 

with again an average lattice spacing oq defined as before. 

But there is a better way to derive correctly the modified form of the wave function. From 
the asymptotic solution for the wave function of Eq. (jll5p it is possible to first obtain a partial 
differential equation for ■0(^tot5 Vtot)- The equation reads (in the following we shall write Rtot as R 
and Vtot as V to avoid unnecessary clutter) 

d'^ip dtp dip d"^ ip , , / ,N 

^ + + + cvugygj^ + cnnjj^, + c,^ + c,ur.^P = . (124) 



It is tempting to try to extract a critical exponent from the result of Eq. p2ip . In analogy to the wave 
functional for a free scalar field with mass m, and thus correlation length ^ = 1/m, one would obtain for the 
correlation length exponent f (with v defined by ^ ~ Is ~ ffcl"") from the above wave function the semi-classical 
estimate 1/ — | . In the 2 + e perturbative expansion for pure gravity one finds in the vicinity of the UV fixed point 
= (d- 2) + |(d- 2)^ + C'((d- 2)^) [33l [M [35] . The above lowest order lattice result would then agree only with 
the leading, semi-classical term. 
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The coefficients in the above equation are given by 



n + 9g N3 ^ n + 9go 22 + 9^0 jVa^' R 2 

'''' 2g2 'v 2ql " F ^ 48 ^/2 3i/3 ^ " ^4/3 + J 

2 ll + 9go iVsi? , ...„2^ 



2^ 
9 



32 A _ 32 4^/2A ii . 

16 i? 16 ii 2 



- + . (125) 



Note that in the small curvature, large volume limit [this is the limit in which, after all, Eq. (IllSh 
was derived] one can safely set the coefficients cr and crr to zero. It is then easy to check that 
the solution in Eq. (jllSp satisfies Eqs. (jl24p and (|125p . up to terms of order 1/^2. Also note that 
here, and in Eqs. (jll2p . (|113p . (|114p and (jllSp . we take the large volume limit F — )• oo, treating 
the number of tetrahedra as a large, fixed parameter. A differential equation in the variable V 
only can be derived as well (with coefficients that are functions of R), but then one finds that the 
required coefficients are not real, which makes this approach less appealing. 
In the hmit R^ Q Eq. (I124p reduces to 

d^i, ll + 9go 7V3 d^l^ ^ 32A 

7^ + o 2 ■ 17 ■ in? + ^ = > (126) 

which is essentially Eq. (j84p in the same limit, with solution given previously in Eq. (j9ip . 



8 Nature of the Wave Function Solution i/j 

In this section we discuss some basic physical properties that can be extracted from the wave 
function solution ip{V, R). So far we have not been able to find a general solution to the fundamental 
Eq. (jl24p . but one might suspect that the solution is still close to a Bessel or hyper geometric 
function, possibly with arguments "shifted" according to Eqs. (jl22p and (jl23p . as was the case in 
2 + 1 dimensions. As a consequence, some physically motivated approximations will be necessary 
in the following discussion. Let us discuss here in detail one possible approach. If one sets the 
troublesome coefficient cvr = in Eq. (jl24p . and keeps only the leading term in cy, then the 
relevant differential equation becomes 

d'^ip dip 

Qy2 + CV -Qy + CX^P + CcUTV W = ^ , (127) 
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with coefficients given in Eq. (jl25p . except that from now on only the leading term in cy and Cx 
will be retained (otherwise it seems again difficult to find an exact solution). Note that the above 
equation still contains an excplicit curvature term proportional to R, from Ccurv Now a complete 
solution can be found in terms of the confluent hypergeometric function of the first kind, i-Fi(a, b, z) 
[36 1 137 1 [38], Up to an overall wave function normalization constant, it is 

P ( (11+990) A^a , W2R_ 



il^iV, R) ~ e 90 G 




(ll+9i?o)Af3 _|_ iV2R 
4^0 9oGVA 



(ll + 9go)A^3 iV2R (ll + 9go)A^3 SiV2\V 



. (128) 



Here again is just a number, given previously in Eq. (|51|) . and A''3 the total number of tetrahedra 
for a given triangulation of the manifold. Note that this last solution still retains three key properies: 
it is a function of geometric invariants (y, R) only; it is regular at the origin in the variable V (the 
irregular solution is discarded due to the normalizability constraint); and finally it agrees, as it 
should, with the zero curvature solution of Eqs. ()9ip and (j92p in the limit i? = 0. 

The above wave function exhibits some intriguing similarities with the exact wave function 
solution found in 2 + 1 dimensions; the difference is that the total curvature R here plays the role 
of the Euler characteristic x there. Let us be more specific, and discuss each argument separately. 
For the arguments of the confluent hypergeometric function of the first kind, i-Fi(a, 6, z), one finds 
again b = 2a for R = 0, with both a and b proportional to the total number of lattice sites, as in 
2 + 1 dimensions [2]. Specifically, here one has 



Re{a) = ^^ + 9go _ 0.5464iV3, (129) 
4^0 



whereas in 2 + 1 dimensions the analogous result is 

Re{a) = IN2. (130) 

The curvature contribution in both cases then appears as an additional contribution to the first 
argument (a), and is purely imaginary. Here one has 

Im{a) = / d^x^R, (131) 



qoVXG 

whereas in 2 + 1 dimensions the corresponding result is 



Im{a) = - [ <fx^R. (132) 

2 V 2 A Gr J 
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Finally, here again the third argument z is purely imaginary and simply proportional to the total 
volume. From the above solution 



8V2A 



whereas in 2 + 1 dimensions 



2V2A 

i 



J<fx^, (133) 
j(fx^. (134) 



G 

Nevertheless we find also some important additional differences with the 2 + 1 result, most notably 
the various gamma-function coefficients involving the curvature i?, which are entirely absent in the 
lower dimensional case, as well as the fact that the critical (UV fixed) point is located at some 
finite Gc here [see Eq. ()119p ]. whereas it is exactly at Gc = in 2 + 1 dimensions [2]. 

Let us now continue here with a discussion of the main properties of the wave function in 
Eq. ()128p . First let us introduce some additional notational simplification. By using the coupling 
g [see Sec. dH and Eq. ()36p ] one can make the above expression for tjj slightly more transparent 

p ( (ll+9(?o)Af3 _L 2iR\ 

^p(V, R) ~ e "OS . — ^ — 2 ^ 

^ ^ rfl (ll+9<?o)Af3 I 2iR\ 

f {n + 9qo)N3 2iR (ll + 9go)jV3 ^iV\ 
4^0 Qog Qog 



We remind the reader that, by virtue of Eq. (I5ip . in all the above expressions qo is just a numerical 
constant, qo = 2tt/ cos~^(|) = 5.1043. Note that for weak coupling the curvature terms become 
more important due to the 1/g^ coefficient. The resulting probability distribution |^/^(1/, i?)|^ is 
shown, for some illustrative cases, in Figures 3,4 and 5. 

One important proviso should be be stated here first. We recall that having obtained an (exact 
or approximate) expression for the wave function does not lead immediately to a complete solution 
of the problem. This should be evident, for example, from the general expression for the average 
of a generic quantum operator 0{g) 

(n(n)\ = ^^1^1^) - Id^^[g]■0{g^,)■mg^J]\' 

^ ^'^^ ^ (^1^) " Id^.[g]-mg.,]\' ' ^ ^ 

where dfj,[g] is the appropriate (DeWitt) functional measure over the three- metric gij. Because of 
the general coordinate invariance of the state functional, the inner products shown above clearly 
contain an infinte redundancy due to the geometrical indinstinguishability of 3-metrics which differ 
only by a coordinate transformation [7]. Nevertheless this divergence is of no essence here, since it 
cancels out between the numerator and the denominator. 
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On the lattice the above average translates into 

- ~ Tmpijww ' ^ ' 

where dfi[P] is the Regge- Wheeler lattice transcription of the DeWitt functional measure [7] in 
terms of edge length variables, here denoted collectively by The latter includes an integration 
over all squared edge lengths, constrained by the triangle inequalities and their higher dimensional 
analogs [30]. Again, because of the continuous local diffeomorphism invariance of the lattice theory, 
the individual inner products shown above will contain an infinte redundancy due to the geometrical 
indinstinguishability of 3-metrics which differ only by a lattice coordinate transformation. And, 
again, this divergence will be of no essence here, as it is expected to cancel between numerator and 
denominator [22j. 

It seems clear then that, in general, the full functional measure cannot be decomposed into a 
simple product of integrations over V and R. It follows that the averages listed above are in general 
still highly non-trivial to evaluate. In fact, quantum averages can be written again quite generally 
in terms of an effective (Euclidean) three-dimensional action 

{^\d{gm = M j d^l[g] d{g,,) exp{-5e//b]} , (138) 

with Seff[g\ = — In and M a normalization constant. The operator 0{g) itself can be local, 

or nonlocal as in the case of correlations such as the gravitational Wilson loop [31]. Note that the 
statistical weights have zeros corresponding to the nodes of the wave function so that S^ff is 
infinite there. |^ 

Nevertheless it will make sense here to consider a semi- classical expansion for the 3 + 1- 
dimensional theory, where one simply focuses on the clearly identifiable stationary points (maxima) 
of the probability distribution obtained by squaring the solution in Eqs. (jl28p or ()135p . In the 
following we will therefore focus entirely on the properties of the probability distribution |V'(V^, i?)p 
obtained from Eq. (jl28p or (|135p . For illustrative purpose, the reader is referred to Figures 3,4 and 
5 below. 

As discussed previously, the asymptotic expansion for the wave function at large volumes implies 
the existence of a phase transition at some G = Gc [see for example Eqs. (I118p and (I119p ]. In 
addition, the explicit solution in Eq. ()135p allows a more precise non-perturbative characterization 



In practical terms, the averages in Eqs. (|136p and (|137|) are difficult to evaluate analytically, even once the 
complete wave function is known explicitly, due to the non-trivial nature of the gravitational functional measure; in 
the most general case these averages will have to be evaluated numerically. The presence of infinitely many zeros in 
the statistical weights complicates this issue considerably, from a numerical point of view. 
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of the two phases. In view of the non-trivial and generally complex arguments of both the gamma 
function and the confluent hyper geometric function, the analytic properties of the wave function, 
and therefore of the probability distribution, are quite rich in features, at least for the more general 
and physically relevant case of non-zero curvature. 

One first notes that for strong enough coupling g the distribution in curvature is fairly flat 
around R = 0, giving rise to large fluctuation in the latter (see Figure 3). On the other hand, for 
weak enough coupling g the probability distribution in curvature is such that values around R = 
are almost excluded, since they are associated with a very small probability. Furthermore, unless 
the volume V is very small, the probability distribution is also generally markedly larger towards 
positive curvatures (see Figure 4). 

In order to explore specifically the curvature {R) dependence of the probability distribution, it 
would be desirable to factor out or remove the dependence of the wave function ip{V, R) on the total 
volume V. To achieve this, one can employ a mean-field- type prescription, and replace the total 
volume V by its average {V). After all, the probability distribution in the volume is well behaved 
at large G [see Sec. ([6])], and does not exhibit any marked change in behavior for intermediate 
G [as can be inferred, for example, from the asymptotic form of the wave function in Eq. (|115|) ]. 
Consequently we will now make the replacement in ^{V, R) 

V {V) = iVg ( K ) = 0.2643 = 0.3738 g , (139) 

V A 

obtained by inserting the result of Eq. (|64p. This replacement then makes it possible to plot the 
wave function of Eq. (I135p squared as a function of the coupling g and the total curvature R only 
(in the following we use again = 10 for illustrative purposes); see Figure 5. One then notes 
that for strong enough gravitational coupling g = \fG the probability distribution is again fairly 
flat around i? = 0, giving rise to large fluctuations in the curvature. On the other hand, for weak 
enough coupling g one observes that curvatures close to zero have near vanishing probability. The 
distributions shown suggest therefore a clearly pathological ground state for weak enough coupling 
9 < Sc [oi G < Gc, see Eq. (|119p ]. with no sensible four-dimensional continuum limit. 

At this point some preliminary conclusions, based on the behavior of the wave function discussed 
previously in Sec. ([7]) and the shape shown in Figures 3,4 and 5, are as follows. For large enough 
G > Gc, but nevertheless close to the critical point, the flatness in the curvature probability 
distribution implies that different curvature scales are all equally important. The corresponding 
gravitational correlation length is finite in this region as long as G > Gc, and expected to diverge 
at the critical point, thus presumably signaling the presence of a massless excitation at Gc [see 
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the argument after Eq. (jl2ip ]. On the other hand for weak enough couphng, G < Gc we observe 
that the probability distribution appears to change dramatically. The main evidence for this is the 
shape of the approximate wave function given in Eq. (|128p . which points to a vanishing relative 
probability for metric field configurations for which the curvature is small i? ~ 0. This would 
suggest that the weak coupling phase, for which G < Gc, has no continuum limit in terms of 
manifolds that appear smooth, at least at large scales. The geometric character of the manifold 
is then inevitably dominated by non-universal short-distance lattice artifacts; no sensible scaling 
limit exists in this phase. 

If this is indeed the case, then the results obtained in the present, Lorentzian, 3 + 1 theory 
generally agree with what is found in the Euclidean case, where the weak coupling phase was found 
to be pathological as well \20\ [21] (it bears more resemblance to a branched polymer, and has 
thus no sensible interpretation in terms of smooth four-dimensional manifolds). In either case, 
the only physically acceptable phase, leading to smooth manifolds at large distances, seems to be 
the one with G > Gc- It is a simple consequence of renormalization group arguments that in this 
phase the gravitational coupling at large distances can only flow to larger values, implying therefore 
gravitational anti-screening as the only physically possible outcome. 
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Figure 3: Wave function of Eq. (jl35|) squared, \tp{V,R)\'^, plotted as a function of the total volume 
V and the total curvature R, for coupling g = \fG = 1 and A^s = 10. One notes that for strong 
enough coupling g the distribution in curvatures is fairly flat around i? = 0, giving rise to large 
fluctuations in the curvature. These become more pronounced as one approaches the critical point 
at gc- 

9 Summary and Conclusions 

In this work we have discussed the nature of gravitational wave functions that arise as solutions of 
the lattice Wheeler-DeWitt equation for finite simplicial lattices. The main results of the paper were 
given in Eqs. (jl24p . (jl28p and (jl35p . While there are many aspects of this problem that still remain 
open and unexplored, we have nevertheless shown that the very structure of the wave function is 
such that it allows one to draw a number of useful and perhaps physically relevant conclusions 
about ground state properties of pure quantum gravity in 3 + 1 dimensions. These include the 
observation that the theory exhibits a phase transition at some critical value of Newton's constant 
Gc [given in Eq. (fTTUD ]. 

The structure of the wave function further suggests that the weak coupling phase, for which the 
coupling G < is pathological and cannot be interpreted in terms of smooth manifolds at any 
distance scale. In view of these results it is therefore not entirely surprising that calculations that 
rely on the weak field, semiclassical or small G expansion run into serious trouble and uncontrollable 
divergences very early on. Such an expansion does not seem to exist if the non-perturbative lattice 
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Figure 4: Same wave function of Eq. (I135p squared, ItpCV, R)\'^ , plotted as a function of the total 
volume V and the total curvature R, but now for weaker coupling g = \fG = 0.5, and still N3 = 10. 
For weak enough coupling g the distribution in curvature is such that values around R = are 
almost completely excluded, as these are associated with a very small probability. Note that, unless 
the total volume V is very small, the probability distribution is markedly larger towards positive 
curvatures. 
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Figure 5: Curvature distribution in i? as a function of the coupling g = \fG. The strong couphng 
relationship between the average volume and the coupling g [Eq. (|64p ] allows one to plot the wave 
function of Eq. (|135p squared as a function of the coupling g and the total curvature R only (we 
use again here A's = 10 for illustrative purposes). Then, for strong enough coupling g = ^/G, the 
probability distribution \ip\'^ is again fairly flat around -R = 0, giving rise to large fluctuations in the 
curvature. The latter are interpreted here as signaling the presence of a massless particle. On the 
other hand, for weak enough coupling g one notices that curvatures close to i? = have essentially 
vanishing probability. The distribution shown here points therefore toward a pathological ground 
state for weak enough coupling g < gc [given in Eq. (jll9p ]. with no sensible continuum limit. 
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results presented here are taken seriously. The correct physical vacuum apparently cannot in any 
way be obtained as a small perturbation of flat, or near-flat, spacetime. On the other hand the 
strong coupling phase does not exhibit any such pathology, and is therefore a good candidate for a 
physically acceptable ground state for pure quantum gravity. It is a simple consequence of standard 
renormalization group arguments that in this phase Newton's constant grows with distance, and 
thus this phase exhibits gravitational anti-screening. 

In the Euclidean lattice theory of gravity in four dimensions it was also found early on \20\ [2T] 
[see |32j for more recent numerical investigations of 4d lattice gravity, including the determination 
of the critical point and scaling exponents] that the weak coupling (or gravitational screening) phase 
is pathological with no sensible continuum limit, corresponding to a degenerate lower dimensional 
branched polymer. The calculations presented here can be regarded, therefore, as consistent with 
the conclusions reached earlier from the Euclidean lattice framework. No new surprises have arisen 
so far when considering the Lorentzian 3 -|- 1 theory, using an entirely different set of tools. 

It is also worthwhile at this point to compare with other attempts at determining the phase 
structure of quantum gravity in four dimensions. Besides the Regge lattice approach, there have 
been other attempts at searching for a non-trivial RG UV fixed point in four dimensions using 
continuum methods. In one popular field theoretic approach one develops a perturbative diagram- 
matic 2 + e continuum expansion using the background field method to two loop order |33i [3^ [35] . 
This then leads to a non-trivial UV fixed point Gc = 0(e) close to two dimensions. Two phases 
emerge, one implying again gravitational screening, and the other anti-screening. In the truncated 
renormalization group calculations for gravity in four dimensions |39[ 140] , where one retains the 
cosmological and Einstein-Hilbert terms, and possibly later some higher derivative terms, one also 
finds evidence for a non-trivial UV fixed point scenario. As in the case of gauge theories, both of 
these methods are ultimately based on renormalization group flows and the weak field expansion, 
and are thefore unable to characterize the non-perturbative features of either one of the two ground 
states. Indeed, within the framework of the weak field expansion inherent in these methods, only 
the weak coupling phase has a chance to start with. It is nevertheless encouraging that such widely 
different methods tend to point in the same direction, namely a non-trivial phase structure for 
gravity in four dimensions. 

Let us add here a few more comments, aimed at placing the present work in the context of 
previous similar calculations for the same theory, namely the treatment of quantum gravity via 
the continuum Wheeler-DeWitt equation. A number of attempts have been made over the years 
to obtain an estimate for the gravitational wave functional ^[g] in the absence of sources. These 
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generally have relied on the weak field expansion in the continuum, as originally done in [10\ 112] . 
Thus, for example, one finds in 3 + 1 dimensions 

nh^^] = N exp|-i j d^\^khj^{\^) /if,^*(k)| , (140) 

where /i^"^(k) is the Fourier amplitude of transverse-traceless modes for the linearized gravitational 
field in four dimensions. The above wave functional describes a collection of harmonic oscillator 
wave functions, one for each of the infinitely many physical modes of the linearized gravitational 
field. 

The above wave functional describes a collection of harmonic oscillator wave functions, one for 
each of the infinitely many physical modes of the linearized gravitational field. As in the case of 
the electromagnetic field, the ground state wave functional can be expressed equivalently in terms 
of first derivatives of the field potentials (the corresponding B field for gravity), without having to 
mention Fourier amplitudes, as 

nh^^] = M exp j d^^j dh ^-^^^^^^^^^ ■ (141) 

It is understood that these expressions represent only the leading term in an expansion involving 
infinitely many terms in the metric fluctuation hij. In an expansion about flat space, the cosmolog- 
ical constant term does not appear either. Since Eq. (|140p is just the leading term in the weak field 
expansion, no issue of perturbative renormalizability appears to this order. Nevertheless, higher 
orders are expected to bring in rather serious ultraviolet divergences which cannot be reabsorbed 
into a simple redefinition of the fundamental couplings G and A. Then the results presented in 
this paper [namely Eqs. (jl24p . (|128p and (|135p ] can be viewed therefore as an attempt to extend 
non- perturb at ively the result of Eq. ()140|) beyond the inherent limitations of the weak field limit. 
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